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ABSTRACT 

We show that coronal magnetic dissipation in thin active sheets that sandwich stan- 
dard thin accretion disks in active galactic nuclei may account for canonical elec- 
tron temperatures of a few x 10 9 K if protons acquire most of the dissipated energy. 
Coulomb collisions transfer energy from the ions to the electrons, which subsequently 
cool rapidly by inverse-Compton scattering. In equilibrium, the proton energy density 
likely exceeds that of the magnetic held and both well exceed the electron and pho- 
ton energy densities. The Coulomb energy transfer from protons to electrons is slow 
enough to maintain a high proton temperature, but fast enough to explain observed 
rapid X-ray variabilities in Seyferts. The ~ 10 9 K electron temperature is insensitive 
to the proton temperature when the latter is > 10 12 K. 
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1 INTRODUCTION 

UV and X-ray observations of active galactic nuclei (AGN) , 
particularly in Seyfert 1 galaxies, indicate that the gravita- 
tional binding energy of a massive black-hole is dissipated 
partly in a cold accretion disk and partly in a hot corona 
above it. Thermal Comptonization of soft UV-radiation in 
the corona leads to the production of a hard X-ray contin- 
uum, some of which is reprocessed by the cold disk (Haardt 
& Maraschi 1991, 1993). In order to explain the different 
ratios of X-ray and UV luminosities in different objects, it 
has been suggested that the corona consists of localized ac- 
tive regions. It is likely that these are produced by mag- 
netic fields in the disk amplified through differential rota- 
tion. When the disk magnetic field builds up significantly, 
buoyancy forces the field out and above the disk, giving rise 
to active regions of high magnetic field. Magnetic field buoy- 
ancy can maintain both magnetically and thermally domi- 
nated disks (Galeev, Rosner & Vaiana 1979, Shibata et al. 
1990; Field & Rogers 1993). 

Magnetic dissipation (for example in a reconnection 
site) may more effectively energize protons than the elec- 
tron. Since it is the electrons which radiate, it is therefore 
plausible that the corona has two kinetic temperatures, that 
of the protons, T p and that of the electrons, T c . These tem- 
peratures depend crucially on the energy exchange rate from 
protons to electrons. 

In this Letter we show that an optically thin corona 
above a thin disk can generate the appropriate electron tem- 
perature to explain the observations of the hard continuum 
of AGN if electrons and protons are coupled via two-body 
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Coulomb collisions. Two-temperature models are usually 
considered in the context of low density ion supported tori 
and advection dominated thick disks (e.g. Shapiro, Light- 
man and Eardley 1976; Rees et al., 1982; Narayan and Yi 
1995). In our context however, although the thin magnet- 
ically active coronal sheets are much less dense than their 
underlying thin dense disk, they can still have a much higher 
density than standard advection dominated disks. We show 
that this enables enough energy transfer by Coulomb elec- 
tron heating to satisfy constraints of rapid X-ray variability 
observed in many Seyfert I nuclei (e.g. Zdziarski et al 1994, 
1995). Also, the active region need not be infalling, but can 
in fact be wind or ejection dominated and still transfer the 
required energy to electrons. 

To obtain this solution, we balance the energy transfer 
rates from the magnetic field to the protons and from the 
electrons to radiation. Balancing these rates does not pre- 
sume that any of the component energy densities are equal. 
We find that the proton and magnetic energy densities are 
both much larger than equilibrium electron and photon en- 
ergy densities. 



2 MAGNETIC ENERGY DISSIPATION 

Magnetic dissipation sites (e.g. reconnection and/or shocks) 
in astrophysical plasmas always involve large gradients in 
the magnetic field, and the presence of length scales much 
smaller than the field gradients of the non-dissipating re- 
gions. This is because the induction equation governing the 
magnetic field 
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dB/dt = V x (v x B) + v m X7 2 B, (1) 

where v m is the magnetic diffusivity, implies flux freezing 
and no dissipation as long as the last term is small. The 
ratio of this last term to the second is the magnetic Reynolds 
number R m = vL/r) » 1, where L and v are the magnitude 
of velocity and gradient length scales. The large sizes and 
conductivities in AGN (and astrophysical systems generally) 
make dissipation negligible except in regions of large field 
gradients, ie. thin current sheets (e.g. Parker 1979). 

The more pronounced the variability of the AGN source, 
the more likely the fewer the number of coronal dissipa- 
tion regions: if N dissipation regions operate at a time, then 
a fluctuation changing the number of dissipation regions, 
would change the luminosity ~ (N ± n) /N, where n is the 
change in number of regions. If N » 1, n ~ N 1/2 so that 
the effective change in luminosity would be small compared 
to when N ~ 1. 

A natural site of coronal dissipation by reconnection 
and its slow shocks, would be at the interface between the 
newly rising field from the disk and any pre-existing coro- 
nal field (Begelman 1990; Shibata et al. 1990). Disspation 
between pre-existing loops has also been considered (Field 
& Rogers 1993). In general, the dissipation layer could be 
expected to cover a sizable cross sectional fraction of the 
underlying disk. 

The development of small turbulent length scales in the 
magnetically active regions enables the transfer of field en- 
ergy into particles (e.g. Larosa et al. 1996). The turbulence 
generates a cascade of magnetic scales and the true dis- 
sipation scale is that below which the energy drains into 
particles rather than into smaller magnetic structures. The 
transfer of magnetic to particle energy likely occurs from 
MHD wave-particle interactions, for example by stochastic 
or shock Fermi energization (e.g. Achterberg 1987; Eilek & 
Hughes 1991; Larosa et al., 1996). Such energization requires 
two types of particle scattering, both of which likely favor 
protons for AGN coronae as we now discuss. 

Particle scattering by MHD waves is the first required 
mechanism. The MHD waves are essentially magnetic com- 
pressions on scales much larger than the particle gyro-radii 
which can "mirror" particles as they approach. For stochas- 
tic Fermi energization, the scatterers move in random direc- 
tions at ~ va, the Alfven speed. Only the "head-on" colli- 
sions between particles and scatterers transfer energy to the 
particles. For particles moving faster than or of order va, 
both head-on and trailing scatterings occur. Even though 
much slower particles do see only head-on collisions, there 
is a critical velocity perpendicular to the magnetic field re- 
quired by particle in order for them to be scattered (Larosa 
et al. 1996). This latter condition dominates, and more en- 
ergy is gained for faster particles. If the protons have higher 
velocities to begin with, then proton rather than electron 
scattering would be favored. 

The above scattering can only enhance a particle's mo- 
mentum component along the field. In order to be repeatedly 
scattered and gain significant energy, the particles must be 
pitch angle isotropized with respect to the local magnetic 
field. This pitch angle scattering is accomplished by reso- 
nance with plasma Alfven waves on the scale of the particle 
gyro-radii (Achterberg 1987; Eilek & Hughes 1991). Such 
waves only exist below a critical frequency and so the par- 



ticles must have significantly large gyro-radii before gyro- 
resonance can occur. The required pitch angle scattering 
also favors the protons: in a plasma for which the Alfven 
speed is of order c, the electrons need to be quite relativis- 
tic to resonate with Alfven waves. Such favoring of protons 
has been called the "electron injection" problem (Eilek & 
Hughes 1991). 

Shock Fermi energization may also operate across slow 
shocks in the reconnection region (Blackman & Field 1994). 
This produces only head- on collisions and is thus more effi- 
cient than stochastic processes. However, similar arguments 
to those above would apply, and protons would be more eas- 
ily energized. 

The previous discussion suggests that when the pro- 
tons are already at higher kinetic temperatures than the 
electrons, magnetic dissipation can maintain such a state in 
equilibrium. Setting up the initial configuration is a separate 
problem since the initial plasma might not necessarily have 
protons at a higher temperature than electrons. However, 
even if the initial dissipation went into electrons, they would 
cool very fast. The proton energy density could slowly grow. 
After some time, the protons could acquire enough energy 
to dominate the electrons, leading to the two-temperature 
equilibrium we describe herein. We leave a more detailed 
investigation for future work. 

Simulations (Priest & Forbes, 1986) as well as solar 
(Tsuncta 1994) and the geomagnetic tail observations (Coro- 
niti et al. 1994) indicate that magnetic reconnection oc- 
curs in the presence of slow shocks. This suggests that 
Petschek (Petschek 1962) type reconnection can occur in 
nature, distinctly from Sweet-Parker (Sweet 1958; Parker 
1979) type reconnection. The main difference between the 
two for non(or mildly)-relativistic Alfven speeds is that for 
the latter, the ratio of height to width of the the overall dis- 
sipation region satisfies is h/r ~ v&\ a /VA ~ 1/Rm 2 while in 
the former the ratio is h/r ~ v^/Va ~ l/log c [R m ] where 
the R m >> 1 is measured outside of the dissipation region. 
The Sweet- Parker dissipation region is generally thinner and 
magnetic energy conversion occurs much slower. Our cal- 
culations below, combined with observational constraints, 
seem to suggest that h/r ~ 1/10 is reasonable, possibly sug- 
gesting that the dissipation regions are Petschek- like. 

Finally, note that the processes discussed above 
might generate non-thermal proton distributions. However, 
the electron energization times are also of order their 
Maxwellian relaxation times (Spitzer 1956). Thus the elec- 
trons can be rapidly thermalized even if the protons arc 
initially non-thermal. 

3 A TWO-TEMPERATURE CORONA 
3.1 Time scales 

We now motivate the model by exploring relevant time scales 
using simple approximate formulae, reserving a fuller treat- 
ment for the next Section. 

Thermal Compton scattering gives the power-law shape 
to the observed hard X-ray spectrum and exponential cut- 
off of Seyfert 1 nuclei, provided that the plasma is optically 
thin to electron scattering (i.e. r < 0.4) and has an electron 
temperature kT c ~ 2 — 5 x 10 9 K. The Compton cooling time 
scale 
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(2) 



where the compactness parameter I — L/r{oT /m c c 4 ). L is 
the luminosity of the active region which is of size r. For 
a source such as MCG-6-30-15 which has a luminosity of 
10 43 ergs _1 and varies typically on a time scale of 1000 s 
then £ ~ 10, so t coo \ ~ icross/3, where t croS s = r/c. The 
electrons cool on less than a light crossing time. 

The electron-proton coupling time, i ep is defined from 
the energy transfer equation (Spitzer 1956) 

dTe (T p - To) 

dt tcp 

Wm p (0c + P ) 3/2 100 



2 m c n <TTclnA tt 

where we have assumed that the sum of e = kT e /m e c 2 and 
P = kT p /m p c 2 is of order unity. Provided that T c <C T p , 
the coupling time to attain T c 



te — te 



0.1 



(4) 



Thus if T c ~ 10 9 K and T p ~ 10 12 K, t e w t cross /4, if 
tt = 0.4, which means that t c ~ i C ooi- This last condition is 
required for heating to balance cooling for the electrons and 
determines T e . 

A slab geometry enables T c ~ 10 9 K for lower values 
of T p . If, for instance, the thickness h is one tenth of the 
width r, the cooling time does not change, but the higher gas 
density means that the coupling time decreases by the factor 
h/r. Moreover, if the slab is favourably oriented so that it 
is perpendicular to the line of sight, a luminosity variation 
can be observed on a time scale of h/c, shorter than the one 
earlier assumed by the same factor. This enables occasional 
variation to be detected on 100 sec (see e.g. Reynolds et al 
1995). Thus T p « 30(h/r)(e/rr)T e . 



3.2 Magnetic Field and Particle Energy Densities 

We take r to be the Thomson optical depth given by r = 
n e (TTh. The size h of an active region in the corona can 
be constrained by observations of the shortest luminosity 
variability time scales h/c w 100 s implying that h ~ 3 x 
10 12 cm. For r ~ 0.3, n ~ 1.5 x lO^cm" 3 . 

Given that a fraction ij of the accretion power is dis- 
sipated in the corona ((1 — 17) in the underlying disk), we 
have 



r\L = 4irr ce 7 



de„ 



v dt ' 



(5) 



where e 7 and e m are the photon and magnetic energies re- 
spectively and r is the radius of the effective dissipation 
sheet of thickness h. Now, 47rr 2 fr(cfe m /Gfe)_ ~ 47rr 2 t>di s x 
B 2 /8n so from eqn. 4 we have 



B /87r ~ e 7 (c/u d is). 



(6) 



We assume that the dissipation velocity v diB is a fraction 
(/ ~ 1) °f the Alfven speed va, i.e. Vdis = fB/^/Airp. In a 
reconnection site though, the fraction / ~ h/r. From eqn. 6 
we then have 



B = 7.9 x 10' 



( yL_ 



\ 10 43 erg s- 



1/3 



10 _13 g cm - 



1/6 



( h - v 

V3 x 10 12 cm/ 



1/3 / 



Vl0r 9 



-1/3 



M 



10 7 M Q 



-1/3 



G, (7) 



where r g = GM/c 2 and M is the central (black hole) mass. 
The characteristic Alfven velocity would then be va ~ 0.2c. 

The above estimate uses only the energy transfer rates. 
We now show that the assumption of equipartition between 
the magnetic energy density e m , and proton energy den- 
sity e p = e m would likely underestimate the proton tem- 
perature. In the steady state we know that |(de m /cK)- 1 = 
\(de p /dt) + \ = \(de p /dt)-\, and the additional assumption 
of e m = e p would mean that t ep = e p /\(de p /dt)-\ = 
em/\(dem/dt)-\ = h/vdi B - But using eq (3) with h in place 
of r, this would imply w d i s ~ rc/100 ~ c/300. Since the 
dissipation physics of section 2 suggests 7Jdi s ~ h/r, and 
variability constraints imply h/r ~ 1/10, the assumption 
of equipartition between the magnetic and proton energy 
densities would under estimate the proton temperature by 
a factor ~ 10 compared to what we find in the next section. 

We will see that the reason for this is that t cp can be 
one or two orders of magnitude longer than t c . Therefore, 
before a steady state can ensue, the proton energy density 
must build up to some critical value. The protons act as 
a growing sea of magnetic energy until they can transfer 
energy to electrons at the rate they gain energy from the 
magnetic field. Most of the energy in the active region thus 
resides in the hot protons. The least energy is contained in 
the electrons, whose energy density is less than the photon 
energy density e 7 by the factor £ C ross/icooi- The long t cp also 
means that although a flare may rise on a crossing time, its 
decay time is uncertain. The luminosity might drop grad- 
ually as the proton energy is exhausted or abruptly if the 
region expands. 



4 A DETAILED MODEL 

4.1 Heating of electrons by ions 

We determine T c by taking into account the relevant cool- 
ing processes. We solve (de e /dt)- = (de y /dt) + , taking 
the transfer for Coulomb collisions between populations of 
Maxwellian distributions using the Rutherford scattering 
cross-section (Stepney & Guilbert 1983). We show that typ- 
ically, T ~ few x 10 9 K and that Coulomb interactions are 
fast enough to explain the observed variability. Future work 
should consider how sensitive the result is to the presence 
non-thermal protons, as the energization mechanisms of sec- 
tion 2 would likely provide them. 
We write 



(A 1M\ 3 m ' 

(dEc/dt)- = n c ni(7TC 



2 m D 



(kTj - fcTc) 
K 2 (l/9 e )K 2 {l/9 i 



■In A 



2(0 C + 0if + 1 



*(^)+«- (^) 



(8) 



where InA ~ 20, is the Coulomb logarithm, the K's are the 
modified Bessel functions, and c ,0i are the dimensionless 
electron and ion temperatures previously defined. (Note that 
eqn. 6 is valid over all electron and proton temperatures). 

Electrons in the corona are cooled through inverse 
Compton scattering off the soft photon field produced by 
the accretion disk, through thermal synchrotron radiation 
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Figure 1. The dimcnsionless electron temperature C plotted 
against ion temperature 7]. The dashed line is for a r = 0.1 i.e. a 
lower density and lower magnetic field case. The synchrotron radi- 
ation and its comptonized emission become less important cooling 
components in this case and the electron temperature is slightly 
higher. The two other lines are for r = 0.3. The dashed-dotted 
line is for a M = 10 Mq as opposed to the other lines which 
correspond to typical AGN black hole mass of M = 1O 7 M . The 
electron temperature is relatively independent of the proton tem- 
perature as long as the latter become mildly relativistic. 



and its Comptonization. Hence, (de- 1 /dt) + — g sync h + <?cs + 
gcomp- The dominant term is usually gcomp = 4:8 e caTnU rll A 
where err is the Thomson cross section and £/ ra d = (1 — 
77) L/4ttR J isk . The thermal synchrotron emission rises steeply 
with decreasing frequency, becomes self-absorbed below i/ c 
and gives rise to a blackbody spectrum.. The cooling is 
strongly peaked around v c and can be approximated by 
? S ynch = (2-k /S)m E 6 e u^ / r (Narayan & Yi 1995) where u c 
is obtained by equating the Rayleigh- Jeans emission to the 
synchrotron one (Zdziarski 1985) The Comptonization of the 
synchrotron radiation (qsc) is calculated using Dermer and 
Liang (1991) approximate treatment. 

Fig. 1 shows the electron temperature solution as a func- 
tion of the ion temperature when (de e /dt)- = (de 1 /dt) + is 
solved. The figure illustrates how Coulomb heating of elec- 
trons in the optically thin plasma of the localized regions 
naturally give rise to a T e consistent with X-ray observa- 
tions of high energy spectral cut-offs in AGN. Such temper- 
atures do not depend strongly on T p as long as protons are 
mildly relativistic. The T c is anticorrelated with r; as r in- 
creases, the density p and consequently the magnetic field 
B increase. This implies more synchrotron and Comptonised 
synchrotron cooling to balance an increased Coulomb cou- 
pling efficiency, resulting in a lower 9 e . The temperature 
solution is not a strong function of the black hole mass M 
(see Fig. 1). 

We note that for 0c and r of interest, electron-electron 
scattering through M0ller cross section competes with vari- 
ous loss mechanisms and in particular with Compton losses 
in our case. The question arises of whether electrons can 
thermalize where the electron-electron relaxation time scale 
t ee = 40F6>e /2 / lnA (V cr ) £ tool- Ghiscllini, Guilbert & 
Svensson (1988) though, found that cyclo/synchrotron self- 
absorption as considered in our model, acts as a very effi- 
cient thermalizing mechanism as long as the magnetic energy 
density dominates the radiation energy density, which is a 
condition postulated here. We therefore expect electrons to 
be able to achieve a Maxwellian distribution. 



4.2 Time scale for electron heating 

In the previous section we have shown that electron-ion 
Coulomb coupling can explain the values of T c ~ few x 10 9 K 
inferred from observations. If particles are coupled only by 
two body collisions, the protons must pass energy to elec- 
trons fast enough to explain the rapid variability observed in 
X-ray sources. This requires t CJ = e 7 /(de e /dt) = t c e 1 /e e has 
to be less than the light crossing time h/c, where (de e /dt)- 
is given by eqn. (8). Now since norc = cr/h, we have 



Figure 2. The function ^ plotted against ion temperature T\ As 
indicated in the figure, the solid line is for 8 C = 0.4 and the the 
dashed one for C = 0.6. For T\ > 10 12 the function * becomes 
virtually constant so as to give a ratio IP/t < 1 for typical electron 
temperature and therefore an electron heating time less than or 
comparable to the light crossing time for the active current sheet 
in the corona. 

where f(0 c , 6i) is the function of C , 9\ in eqn. (8) and t e is ex- 
pressed as a fraction ($) of the time scale on which photons 
escape (h/c) a dissipation region. In eqn. (12) ^ < 1/10 
and almost invariant of 7] (for 7- > 10 12 K (sec Fig. 2), 
so that for a typical r < 0.4, t e < h/c, consistent with 
ta-y/t c ~ c/h ~ tcE^/so equal to a few. Variability con- 
straints are thus consistent with ion-electron Coulomb col- 
lisions being the dominant process for electron heating. 



5 DISCUSSION 

We have demonstrated that an electron temperature of a 
few xlO 9 ^, determined from interpretation of X-ray ob- 
servations, can be very naturally explained by an optically 
thin two-temperature plasma corona model with TJ >> T c if 
Coulomb scattering is the dominant source of electron heat- 
ing. If magnetic dissipation in the corona preferentially ener- 
gizes the ions, then T\ will naturally exceed T e . The efficient 
cooling of electrons then leads a T c being locked at the re- 
quired values. The result of T e ~ 10 9 K is robust because T e 
is relatively insensitive to T p when the latter exceeds 10 12 K. 

We have also shown that the time scale for electron 
heating by Coulomb coupling is consistent with the con- 
straints required by observations of fast variability in AGN. 
Though we have considered primairly thermalized ions, the 
actual dissipation may lead to non-thermal ion popula- 
tions. This would not likely effect t e strongly, but the effect 
should be considered in future work. Regardless of the pro- 
ton distribution, the electron distribution would still become 
Maxwellian. 

Pronounced variability suggests that the corona likely 
dissipates in large flares, rather than in a large number of 
small dissipation sites. Since the typically observed variabil- 
ities range between 100 — lOOOsec and not much shorter, 
the effective aspect ratio of the dissipation region can be 
constrained to satisfy h/r ~ 100/1000 = 1/10. We can spec- 
ulate that a small number of operating current sheets would 
mean that this aspect ratio characterizes the dimensions of 
a typical reconnection region itself. Such a thick aspect ratio 
might suggest that Petschek rather than Sweet-Parker type 
reconnection is occurring. In any case, the physics of dissi- 
pation could in principle constrain this value. Future work 
is needed to fully address the dissipation physics. 

The model predicts that the proton energy probably 
dominates the magnetic energy density, but both are much 
larger than the photon and electron energy densities. The 
electron energy density is less than the photon energy den- 
sity by a factor of the cooling time over the light crossing 
time. These differences in energy contents highlight the fact 
that balancing the rates of energy transfer does not neces- 
sarily imply equipartition between any energy densities. 
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Throughout this paper, we have assumed that the only 
coupling between ions and electrons is via Coulomb colli- 
sions which leads to a two-temperature plasma. However, 
Begelman & Chiueh (1989, hereafter BC) have discussed 
a plasma instability which may lead to thermal coupling 
of ions and electrons in a two-temperature flow (71 > T c ) 
on a time scale shorter than the Coulomb coupling time. 
The instability may grow in regions with a high level of 
small scale MHD turbulence that have a plasma (5 parame- 
ter /3 = 8nnkTi/B 2 > 1. We can estimate the rate of electron 
heating via the BC mechanism using the specific parameters 
which characterise our solution and compare it to the rate 
of energy transfer via Coulomb collisions (eqn. 8). For sim- 
plicity we choose 71 = 10 12 . According to BC, the rate of 
heating of electrons when their instability is fully developed 
is given by 



\dt) BC ~ ; A Di \ c ) \L ) 



(10) 



where / is the filling factor of the plasma where the strong 
turbulence needed for the instability is present, mi is the ion 
mass, Wti is the ion thermal velocity, L is the local density 
gradient scale and Aoi is the ion Debye length and A c i its 
respective gyroradius (where A c i = {P/0i)\t>{). We find that 
for va/c ~ 1/5 and for L satisfying the condition for the BC 
instability to operate (eqn. 5.6 in BC paper) with our choice 
of P, we have that 



<0.3/ 



(11) 



(de/dt) B i 
(de/dt) c 

This shows that even in the occurrence of the BC instability, 
the Coulomb heating usually dominates for the choice of 
parameters of interest in our solution. 

Finally, since the coronal plasma discussed herein has 
very high ion temperatures, e.g. 71 > 10 12 we expect col- 
lisions between protons in the high energy tail of the dis- 
tribution to produce pions These would then decay into 
high energy ~ 70 MeV7-rays. Such 7-ray emission would 
constitute a testable feature for this model. Power-law or 
Maxwellian proton (or combination of the two) distributions 
could be distinguished by the different 7-ray spectra they 
predict (Dermer 1986). 
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* This process cools the protons without heating the electrons 
and for ion-temperature of interest it would not cool the ions 
faster than Coulomb collisions. 
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